Detecting the impacts of humidity, rainfall, temperature, and season on chikungunya, dengue and Zika viruses in Aedes albopictus mosquitoes from selected sites in Cebu city, Philippines

Background Aedes albopictus is the secondary vector for dengue virus (DENV) in the Philippines, and also harbors chikungunya (CHIKV) and Zika (ZIKV) viruses. This study aimed to determine the minimum infection rates (MIRs) of CHIKV, DENV serotypes, and ZIKV in Ae. albopictus collected from selected two-site categories by altitude (highland [H] and lowland [L] sites) in Cebu city, Philippines during the wet (WS) and dry seasons (DS) of 2021–2022, and to explore the relationships between these arboviral MIRs and the local weather. Methods The viral RNA extracts in pooled and reared adult Ae. albopictus collected during the DS and WS from two-site categories were subjected to RT-PCR to amplify and detect gene loci specific for CHIKV, DENV-1 to DENV-4, and ZIKV and analyzed with the weather data. Results The range of CHIKV MIRs was higher in the WS (13.61–107.38 infected individuals per 1,000 mosquitoes) than in the DS (13.22–44.12), but was similar between the two-site categories. Rainfall (RF) influenced the CHIKV MIR. The MIR ranges of both DENV-2 (WS: H = 0, L = 0; DS: H = 0–5.92; L = 0–2.6) and DENV-4 (WS: H = 0, L = 0–2.90; DS: H = 2.96–6.13, L = 0–15.63) differed by season but not between the two-site categories. Relative humidity (RH), RF, and temperature did not influence DENVs’ MIRs. The MIR range of ZIKV was similar in both seasons (WS: 11.36–40.27; DS: 0–46.15) and two-site categories (H = 0–90.91, L = 0–55.56). RH and temperature influenced ZIKV MIR. Conclusions RF influenced CHIKV MIR in Ae. albopictus, whereas RH and temperature influenced that of ZIKV. Season influenced the MIRs of CHIKV and DENVs but not in ZIKV. Ae. albopictus were co-infected with CHIKV, DENVs, and ZIKV in both highland and lowland sites in Cebu city. Recommendations include all-year-round implementation of the Philippine Department of Health’s 4S enhanced strategy and installation of water pipelines in rural highlands for vector and disease control. Our findings are relevant to protect public health in the tropics in this climate change. Supplementary Information The online version contains supplementary material available at 10.1186/s12985-024-02310-4.


Background
While Aedes aegypti (Linnaeus) is the main vector of arboviruses such as chikungunya (CHIKV), dengue (DENV), and Zika (ZIKV), Aedes albopictus (Skuse), the Asian tiger mosquito, plays also a role in its prevalent arboviral transmission to humans [1][2][3][4] throughout the tropics and subtropics [5][6][7][8].Dengue fever is the most common global arboviral infection; 3.6 billion people reside in dengue-endemic areas [9].Chikungunya reached its global significance in the 21st century [10], whereas Zika caused 7.6 million cases globally in 2016 [11].In the Philippines, chikungunya infections (600 cases; 0 death) were 545% higher from January 1 to December 31, 2022 than in the same period in 2021.Of which, 127 cases occurred in Central Visayas, where Cebu city belongs [12].The Philippine dengue surveillance (61,170 cases; 216 deaths; case fatality rate [CFR] = 0.4%) ranked second and first in the number of cases and deaths, respectively, in the Western Pacific Region in 2021 [13].Of which, 2,604 cases and 17 deaths (CFR = 0.7%) occurred in Central Visayas during the same period of 2021 [12].The same trend in 2022 was reported in the Philippines (205, 679 cases; 672 deaths; CFR = 0.3%) [14]; of which 17,894 dengue cases (104 deaths; CFR = 0.6%) were reported in Central Visayas [12].Meanwhile, the first 17 Zika cases were reported in the Philippines in 2016 without travel history from an affected country a month prior to their onset of illness [15].No Zika case was reported since 2016, but because Ae. aegypti and Ae.albopictus are found in the country, CHIKV and ZIKV might have been circulating around.
Most studies revealed that temperatures influence large epidemics of mosquito-borne arboviral diseases [16][17][18].In Latin America, temperature and its range, rainfall (RF), and population size are the important predictors of Zika transmission [19].In Brazil, RF predicts ZIKV and CHIKV with a lead time of three weeks each time [20].In Mexico, the dynamics of dengue, chikungunya, and Zika are strongly associated with seasonal climatological variability and transmission potential of these arboviral pathogens [21].In Asia, only a few studies explored the relationship between arboviral mosquito-borne diseases and the weather.In India, Shil et al. [22] observed a strong positive association between the occurrence of CHIKV cases and RF variations.In Japan, Furuya's model [23] predicts that if the daily mean temperature rises from 28 to 30 °C, the median arboviral reproduction number (R 0 ) increases by 18% for dengue, 4.3% for chikungunya, and 11.1% for Zika.However, most studies monitored only CHIKV and ZIKV in Aedes such as in Asia-Pacific [24], Hindu-Kush Himalayan region [25], Singapore [26], Sri Lanka [27], and Thailand [7].In the Philippines, only Balingit et al. [28] and Edillo et al. [29] reported DENVs in Ae. aegypti, the primary dengue mosquito vector, in Tarlac and Cebu cities, Philippines, respectively.However, no study, to our knowledge, has determined the influence of the weather to the infections of CHIKV, DENVs, and ZIKV in the Philippine Ae. albopictus, the secondary dengue mosquito vector, in this current climate change.
The National Dengue Prevention and Control Program of the Philippine Department of Health (DOH) [30] includes the following components: (1) surveillance through the Philippine Integrated Disease Surveillance and Response, (2) case management and diagnosis, (3) integrated vector management, (4) outbreak response, (5) health promotion and advocacy, and (6) research.The enhanced 4 S strategy [12] has been the main focus for dengue prevention and control, where "4S" stands for the following: (1) seek and destroy mosquito-breeding sites, (2) seek early consultation if one develops dengueassociated symptoms, (3) employ self-protection measures such as wearing long pants and long-sleeved shirts, and (4) say "no" to indiscriminate fogging, and implement fogging only during outbreaks in hotspots.Climate change affects the distribution of Ae. aegypti and Ae.albopictus, which allows them to spread new pathogens to initially-unaffected populations in the mountains [31].The latter calls for urgent actions on vector control strategies [32], including source reduction of vectors, to manage the spread of these arboviral diseases, and to protect public health [33][34][35] considering the lack of effective vaccines available.
We hypothesized that weather conditions might influence CHIKV, DENVs, and ZIKV in Ae. albopictus from selected highland and lowland sites in Cebu city, Philippines during the wet (WS) and dry seasons (DS).Thus, this study aimed: (1) to determine the presence and calculate the minimum infection rates (MIRs) of CHIKV, DENV-1 to DENV-4, and ZIKV in Ae. albopictus collected from selected two-site categories by altitude (highland and lowland sites) in Cebu city during the 2021-2022 WS and DS, and (2) to explore the relationships among the arboviral MIRs with the local weather in the study sites.Ae. albopictus was chosen in this study for the following reasons: (1) there has been no study in the Philippines, to our knowledge, that reported the impacts of weather, season, and altitude to the arboviral infections in this mosquito species.(2) Arboviral monitoring in Ae. albopictus is a relevant strategy for tracking their activities and for implementing direct control actions [28,30,36] against associated diseases and vectors as influenced by the weather.However, the local Philippine DOH focuses in monitoring the number of suspected arboviral mosquito-borne diseases in humans and not in mosquito vectors, and community participation of the enhanced 4S strategy.(3) Also, Ae. albopictus is a good sentinel species in DENV prevalence during the interepidemic period [37] in this climate change.

Study sites
The study sites in Cebu city, Philippines were selected barangays (i.e., smallest government units) based on: (1) altitude (highland vs. lowland), (2) presence of local dengue cases, (3) vegetation cover, and (4) in-between site distance of more than 500 m (Table 1; Fig. 1).Highland sites had elevations between 400 and 700 m above sea level (m ASL); lowland sites, between 60 and 80 m ASL.

Mosquito collections
Aedes eggs and sub-adults were collected bimonthly during the DS (March-May 2021 and 2022) and WS (June-November 2021; February and June 2022) using the modified ovicidal/larvicidal (O/L) traps or ovitraps (Department of Science and Technology, Manila).Briefly,  the ovitrap consisted of a black plastic container, filled halfway with water, and a filter paper-wrapped wooden panel.The original O/L pellets were discarded to allow the oviposited Aedes eggs and sub-adults to survive [34,35].Each highland and lowland site had 10 modified O/L traps, placed in locations with thick vegetation cover.The traps were left for 10-15 days.The filter papers with oviposited Aedes eggs (referred to as "egg papers" henceforth) were collected, dried, placed in a clean and dry container and stored in a dark place.Sub-adults were collected from the traps or other water-filled containers within the vicinity by scooping with nets or aspirating with Pasteur pipettes.

Rearing of collected Aedes eggs and sub-adults
Field-collected Aedes egg papers and sub-adults were transported to the University of San Carlos -Mosquito Research Laboratory (USC-MRL), Biology Department, Talamban campus, Cebu city, Philippines.Aedes egg hatching was performed by submerging the egg papers in 10% ascorbic acid solution [38].Aedes larvae and pupae were transferred to a clean plastic container filled halfway with distilled water (DW) and covered with a finemesh cloth.They were fed daily with a dash (~ 0.02 g) of fish food (Sakura; All Aquariums Co., Ltd., Thailand).The DW in the larval plastic container was changed regularly to prevent the formation of bacterial scum until adult emergence.Adult Ae. albopictus samples were identified [39] and transferred to a separate container.They were fed with 10% sucrose solution for five days, individually transferred to a microcentrifuge tube, and stored in a -80 °C ultralow freezer until further processing.

Sample size
The sample size was calculated to consider the difference between binomial proportions (i.e., detection or non-detection) [40] of CHIKV, DENVs, and ZIKV in Ae. albopictus collected from highland and lowland sites of Cebu city during the DS and WS.Neyman allocation was used as an adjunct for generating a sample size for each site (i.e., apportioned from computed total sample size) with an allocated fixed budget [29].At the significance level (α) of 0.05, 80% power (β = 0.2), probabilities p1 = 0.10 and p2 = 0.30 (based from the detection of Guedes et al. [41] and Medeiros et al. [36], respectively) and 20% attrition rate, four pools (i.e., a pool ≤ 30 adults) of Ae. albopictus collected from each site monthly for DS and three pools collected from each site monthly for WS were considered sufficient to warrant rejecting the hypothesis of no difference.A total of 288 pools of Ae. albopictus were prepared for both the DS (n = 144 pools) and the WS (n = 144 pools).

RNA extraction
RNA extraction of adult Ae. albopictus was performed using the RNeasy mini kit (Qiagen, Germany) following the manufacturer's protocol.The quality of RNA extraction was determined by gel electrophoresis.RNA extracts of Ae. albopictus were expected to produce 18s and 28s ribosomal RNA bands at 1,000 and 1,500 base pairs (bp), respectively [42].All extracted sample pools were composed of either male or female, or mixed adults.

Detection of CHIKV, DENV, and ZIKV Reverse transcription-polymerase chain reaction (RT-PCR)
One-step RT-PCR was performed to detect these arboviruses in pooled Ae. albopictus samples using Qiagen kit (Germany; Cat No. 210212).Primer sequences are listed in Additional file 1: Table S1.

Nested PCR of CHIKV
CHIKV nested PCR was done to further amplify the inner region of 6 K/E1 gene.A nested PCR concoction (10 µL) was prepared: 5 µL 2x PCR mastermix, 0.5 µL of both inner 6 K/E1 IA and 6 K/E1 IB primers, 1 µL cDNA template and 3 µL RFW.The PCR profile was set at 95˚C incubation for 10 min, 35 cycles of 95˚C denaturation for 30 s, 55˚C annealing for 30 s, 72˚C extension for 30 s, and 72˚C final extension for 10 min.CHIKV outer and inner primers used in nested RT-PCR were based on Laras et al. [43].

PCR product purification and sequencing of positive controls
Positive controls for CHIKV and ZIKV amplicons were isolated from local Ae.albopictus, whereas those for DENV-1 to DENV-4 were isolated from Ae. aegypti or Ae.albopictus (Table 2).They were subjected to the aforementioned RNA extraction and RT-PCR protocols.For confirmation of arboviral amplicons, they were electrophoresed on 1% agarose gel at 200 VDC for 20 min.
Non-specific amplifications were observed among the PCR products of positive controls, thus, gel excision and extraction of the bands corresponding to expected band size of each arbovirus (Table 2) were done using QIAquick PCR and gel cleanup kit (Qiagen, Germany; Cat.No. 28704) following the manufacturer's protocol.Another PCR run was done to further amplify the desired DNA fragment using Type-it microsatellite PCR kit (Qiagen, Germany).A PCR reaction (20 µL) was made of 10 µL 2x PCR mastermix, 1 µL (10 µM) each of the forward and reverse primers, 4 µL DNA template, and 4 µL RFW.The PCR profile for CHIKV positive control was set at 95˚C incubation for 10 min, 35 cycles of 95˚C denaturation for 30 s, 53˚C annealing for 30 s, 72˚C extension for 30 s, and 72˚C final extension for 10 min.For DENVs and ZIKV positive controls, the PCR profile was similar except annealing temperatures (i.e., 55˚C for DENV and 57˚C for ZIKV).
Purified PCR products for CHIKV, DENV-1 to DENV-4, and ZIKV positive controls were brought to the Philippine Genome Center, University of the Philippines -Diliman, Quezon city for capillary sequencing.Their results were processed using DNA Baser Assembler v4, and were trimmed and cleaned using BioEdit v7.2 [46].The final sequences were subjected to the Basic Local Alignment Search Tool (BLAST) of the National Center for Biotechnology Information (NCBI) [47] for identity verification.

Weather data
The mean monthly relative humidity (RH), monthly total rainfall (RF), and mean monthly temperature nearest (11.8-32.4km) to the study sites throughout the study period were obtained from the local satellite office of the Philippine Atmospheric Geophysical and Astronomical Services Administration (PAG-ASA), Cebu city.

Data and statistical analyses
The monthly MIRs of CHIKV, DENV-1 to DENV-4, and ZIKV were computed for each site by dividing the total number of positive pools by the total number of individual mosquito samples, multiplied by 1,000 [48].Because of non-normality of data distributions, the dependent variables (arboviral MIRs) and the independent variables (season, altitude, and weather data) were subjected to Kruskal-Wallis (K-W) statistical tests with arbitrarily categorized ranks for weather data (Table 3) by using SPSS Advanced Professional (v28).Spearman's rank coefficient  of correlations were performed between the MIRs of CHIKV and ZIKV, and the weather with significant K-W test results to assess the correlation.Pearson's productmoment correlation tests were performed between the MIRs of DENVs and the dengue cases.Two-tailed Student's t-tests were determined for the statistical differences between seasonal weather data.

Isolation of CHIKV, DENVs, and ZIKV in Ae. albopictus as positive controls
The positive controls isolated from the local Ae.albopictus formed the expected gel bands consistently at 200 bp for CHIKV [43], 511 bp for DENV [44], and192 bp for ZIKV [45] (Table 2).The NCBI BLAST results of amplicon nucleotide sequences of these bands yielded 100% identity matches.Likewise, the positive controls with their expected gel bands for DENV-1 (482 bp), DENV-3 (290 bp), and DENV-4 (392 bp) isolated from the local Ae.aegypti, and for DENV-2 (119 bp) isolated from the local Ae.albopictus were formed consistent with Lanciotti et al. [44].5) revealed that seasonal MIRs of DENV-2 (K-W H-test, H = 4.097, df = 1, P = 0.043) and DENV-4 (K-W H-test, H = 11.928,df = 1, P = 0.001) differed, however, all DENV serotypes when grouped in two-site categories by altitude did not differ (K-W H-test, P > 0.05).Results suggest that Ae. albopictus circulated all DENV serotypes in both highland and lowland sites, and their residents were exposed to risks of contracting dengue.Moreover, the seasonality of DENV-2 and DENV-4 in Ae. albopictus across sites may suggest that local weather might have strongly influenced their viral replications.

MIRs of DENVs in Ae. albopictus and dengue cases in Cebu city
While the combined MIRs of all DENV serotypes in Ae. albopictus were high during the DS (March to May 2021 and 2022) and decreased during the start of the WS (July 2021 and June 2022), dengue cases in Cebu city (Additional file 3: Table S2) gradually increased from March 2021 until June 2022 (Fig. 4).Results showed a low, negative, and insignificant correlation (Pearson's correlation coefficient, r (14) = -0.131,P = 0.654), implying a weak association between the MIRs of DENVs in Ae. albopictus and reported dengue cases.Results suggest that the local Ae.albopictus maintained DENVs more frequently during the DS, and transmitted them to humans during the WS.

MIRs of ZIKV in Ae. albopictus by season and altitude
During the DS, the highest MIR (36.81) of ZIKV in Ae. albopictus was recorded in Babag 1 highland site, which was also the highest MIR in the entire sampling period (Table 4; Additional file 2: Dataset 1).The lowest MIR (17.75) of ZIKV was recorded in Babag 2, which was also the lowest MIR in the entire study period.During the WS, the highland Taptap site had the highest MIR (31.95) of ZIKV, while that of Babag 2 was lowest (17.91).
The seasonal mean MIRs of ZIKV in Ae. albopictus did not differ (K-W H-test, H = 1.377, df = 1, P = 0.241) with 45.14 and 38.99 for the WS and the DS, respectively.Likewise, the mean MIRs of ZIKV were similar (K-W H-test, H = 0.000, df = 1, P = 0.993) between the two-site categories by altitude (Table 5).6).By Student's t-test, the seasonal monthly RH (t-test, t (11) = 2.201, P = 0.687) and temperatures (t-test, t (9) = 2.262, P = 0.862) did not differ.Likewise, the seasonal mean monthly total RF did not differ (t-test, t (6) = 2.447, P = 0.739) in the WS and DS of 2021 and 2022.Expectedly, the mean monthly total RF was higher in the WS of 2021 (162.40 mm) than that of the DS of the same year (87.60 mm), however, the mean monthly total RF in the DS (2022) (180.13 mm) was higher than that

Relationship between the MIRs of ZIKV in Ae. albopictus and the weather
Figure 2 shows the monthly ZIKV's MIR plotted against monthly readings of RH, rainfall, and temperature in Cebu city.ZIKV's MIR in Ae. albopictus differed between categorized ranges (Table 3) for RH (K-W H-test, H = 18.106, df = 2, P < 0.0001) (Table 7) and temperatures (K-W H-test, H = 7.067, df = 2, P = 0.029) (Table 7) but not for RF (K-W H-test, H = 1.630, df = 2, P = 0.443) (Table 7).ZIKV's MIR was highest in the high RH category (79-80%) with a mean MIR of 53.42, followed by the low monthly RH category (74-76%) with a mean MIR of 34.22, and lowest at moderate monthly RH (77-78%) with a mean MIR of 29.60.A significant positive Spearman's correlation (Spearman's correlation coefficient, rs = 0.65, P = 0.012) between the MIR of ZIKV and RH was observed.On the other hand, ZIKV's MIR was highest in the low temperature category (27 °C) with a mean MIR of 52.58, followed by moderate (28 °C), and high (29 °C) categories with mean MIRs of 44.27 and 30.47, respectively.A negative and insignificant correlation (Spearman's correlation coefficient, rs = -0.44,P = 0.117) was observed between the MIR of ZIKV and temperature.

Relationship between the MIRs of DENVs in Ae. albopictus and weather
Figure 3 shows the monthly MIRs of DENVs plotted against monthly readings of RH, rainfall, and temperature in Cebu city.K-W test results showed no differences (K-W H-test, P > 0.05) (Table 7) in all DENV serotypes across RH, rainfall, and temperature.

Discussion
This study is the first report in detecting the impacts of Philippine weather (i.e., RH, RF, and temperature) and seasonal (i.e., WS and DS) conditions on the MIRs of CHIKV, four DENV serotypes, and ZIKV in Ae. albopictus, the secondary dengue mosquito vector in the country, collected from selected highlands and lowlands in Cebu city, Philippines during the WS and DS (2021-2022).

CHIKV in Ae. albopictus
First, the seasonal MIRs of CHIKV in Ae. albopictus differed but not in two-site categories by altitude (Table 5).
The former had higher MIR in the WS (2021-2022), suggesting higher CHIKV prevalence during the WS as supported by the significant difference (Table 7) of CHIKV's MIRs in Ae. albopictus as influenced by RF. Results were consistent with reported positive effect of RF on CHIKV transmission [20,[49][50][51].However, similar MIRs of CHIKV in Ae. albopictus in both Cebu city's highland and lowland sites (Table 5) indicate that CHIKV naturally circulated in both two-site categories by altitude.Hence, residents were exposed to risks of contracting CHIKV as evidenced by reported cases described above.Results were consistent in Albanian infestation by CHIKVinfected Ae. albopictus at high altitudes [52].Second, a significant difference of CHIKV's MIR in Ae. albopictus was observed in RF with a weak positive yet insignificant correlation (Spearman's correlation coefficient, rs = 0.16, P = 0.573), but not with the RH and temperature (Table 7).The mean monthly total RF (180.13 mm) during the DS (2022) was higher than that (162.4 mm) of the WS (2022), and huge difference of mean monthly total RF between the DS of 2021 (87.6 mm) and 2022 (180.13 mm) (Table 6).Increased precipitation supports the maintenance of breeding sites contributing to mosquitoes' survival [53].Currently, the correlation between CHIKV's prevalence in Ae. albopictus and RF is consistent with the reports in Thailand [49,51] and Kenya [50].CHIKV is maintained throughout all seasons in Kenya but the abundance of CHIKVinfected mosquitoes is more evident during increased RF (~ 250 mm) [50].Heavy RF with a three-week interval predicts the onset of CHIKV outbreaks [20].Although, this study did not observe significant relationship between CHIKV's MIR in Ae. albopictus and mean monthly temperature nearest to the study sites, Mercier et al. [52] reported that Ae. albopictus only transmits CHIKV at 28 °C.

DENVs in Ae. albopictus
Currently, Ae. albopictus carried all four DENV serotypes (Table 4), like Ae. aegypti, consistent with other studies [54][55][56][57][58], although its odds of having DENV-infected saliva is lower than that of Ae. aegypti [58].Interestingly, seasonal MIRs of DENV-2 and DENV-4 in Ae. albopictus differed (Table 5), suggesting a specific vector-virus interaction.However, the lack of MIR differences of the four DENV serotypes in Ae. albopictus from two-site categories by altitude (Table 5), implies their expanded distribution to highlands as affected by the changing climate [31].Among the highland sites, Taptap had the highest altitude (719 m ASL), where householders stored water for domestic use in artificial containers [34,35].This practice increased the breeding sites in highlands, and may amplify the risks of contracting mosquito-borne diseases.The presence of Ae. albopictus in both Cebu city's lowland (urban) and highlight (rural) sites elucidates their ability to spread these arboviruses [37].With expanded distribution of Ae. albopictus to highlands [59], current vector control programs should be extended there as well.
Currently, the MIRs of DENV-1 and DENV-3 in Ae. albopictus were similar with those of Ae. aegypti [60], whereas those of DENV-2 and DENV-4 were different.In Brazil, DENV-1 and DENV-3 are most common in Aedes and in humans [61].In Singapore, DENV-1 is most common in both Ae.aegypti and Ae.albopictus [62].These studies corroborate with the higher MIRs of DENV-1 and DENV-3 in Ae. albopictus from lowland and highland, respectively (Table 4).
Moreover, Ae. albopictus showed higher MIRs of DENVs (mean = 13.33)during the DS, the inter-epidemic period, than during the WS (mean MIR = 8.85) (Table 4), suggesting the silent DENVs' circulation in mosquitoes on areas with low dengue cases in the DS [63].In addition, DENV-2 was only detected during the DS regardless of altitude.Hence, Ae. albopictus is a good sentinel species in monitoring DENV prevalence during interepidemic period [37].A low negative insignificant correlation (Pearson's correlation coefficient, r (14) = -0.131,P = 0.654) between the MIRs of DENV-infected Ae. albopictus and dengue cases (Fig. 4) could be attributed to the fact that this species is the secondary dengue vector in the Philippines.
Third, the mean MIR of ZIKV in Ae. albopictus did not differ by season and by two-site categories (Table 4).The latter result was consistent between Ae. albopictus populations at highest and lowest altitudes in Albania [52].Current results may suggest that Ae. albopictus transmitted ZIKV all throughout the year and in both lowland and highland sites in Cebu city.Moreover, results imply a whole-year round implementation of the DOH enhanced 4 S strategy in this climate change as in related studies [34,35].

Conclusion
In conclusion, the MIRs of CHIKV in Ae. albopictus were affected by RF, whereas those of ZIKV were affected by RH and temperature.DENVs in Ae. albopictus varied seasonally; DENV-2 and DENV-4 differed seasonally across selected highland and lowland study sites in Cebu city.Seasonal temperature, RF, and RH did not influence the MIRs of DENVs in Ae. albopictus across the two-site categories by altitude.One limitation of this study was the use of PAG-ASA weather data nearest to the study sites but were not measured on-site.However, this study aimed to explore the relationships between the MIRs of the arboviruses in the local study sites with the weather data, hence, an encompassing weather data of Cebu was assumed sufficient to demonstrate such relationships.
This study recommends for an all-year round implementation of the Philippine DOH's enhanced 4 S strategy in this climate change because of the lack of the seasonal impact on ZIKV, DENV-1 and DENV-3 in Ae. albopictus; similar as in our previous studies [34,35].Water pipelines [34,35] and vector control services are highly recommended in the highlands and not just in the lowlands to minimize the practice in storing water in containers that serve as potential breeding sites of dengue mosquitoes as mentioned in the above results.The MIR data of these three arboviruses in Ae. albopictus will be useful in future arboviral modeling at which weather conditions can potentially be used as indicators in predicting possible arboviral outbreaks.Thus, findings of this study are very relevant not just in the Philippines but also in the tropics and subtropics.

Table 1
Coordinates and elevation of the six highland and lowland sites in Cebu city

Table 2
Positive control profiles of the three arboviruses isolated from the local Aedes species

Table 3
Categorized ranks for weather data (mean monthly relative humidity [RH], monthly total rainfall [RF] and mean monthly temperature) utilized in Kruskal-Wallis statistical tests

Table 4 (
sampling period; likewise it had also the lowest MIR (40.00) in the WS.The lowland Pit-os had the highest MIR (57.97) of CHIKV during the WS and in the entire sampling period.K-W test results (Table5) showed that seasonal MIRs of CHIKV in Ae. albopictus differed (K-W H-test, H = 14.814, df = 1, P < 0.0001), with the mean MIR ranks at 30.69 for the DS and at 51.35 for the WS, but did not Additional file 2: Dataset 1) shows that during the DS, the highest MIR of CHIKV (47.45 per 1,000 mosquitoes in the wild; from here onward MIR refers to an estimate for every 1,000 Ae. albopictus in the wild) was recorded in the lowland Talamban site.The lowland Bacayan had the lowest MIR(15.63)in the DS and in the entire

Table 4
MIRs of CHIKV, DENV and ZIKV with number of positive pools in Ae. albopictus by season (dry and wet) and two-site categories by altitude (highland and lowland) in Cebu city, Philippines a n = total number of individuals b number of positive pools out of 24 pools per site

Table 5
H = 0.007, df = 1, P = 0.932) in the two-site categories by altitude.The mean MIR ranks of CHIKV MIR were 42.27 for highland and 42.73 for lowland.Residents in highland sites stored water in artificial containers for domestic use that served as potential breeding sites for Aedes owing to the lack of installed water pipelines by the local government.
Kruskal-Wallis statistical test results for the minimum infection rates (MIR) of chikungunya (CHIKV), dengue (DENV), and Zika (ZIKV) in Ae. albopictus when grouped by season (dry and wet) and by two-site categories by altitude (highland and lowland) in Cebu city, Philippines *Statistically significant differ (K-W H-test,

Table 4
), Ae. albopictus in both two-site categories had the four DENV serotypes.Results suggest that Ae. albopictus played an important role in DENV maintenance during the DS when breeding sites were scarce.K-W test results (Table

Table 6
Weather data (mean monthly relative humidity [RH] and temperature, and monthly total rainfall [RF]) in Cebu city, Philippines during the dry (March -May 2021; March-May 2022) and wet (June-November 2021; February and June 2022) seasons

Table 7
Kruskal-Wallis statistical test results between the minimum infection rates (MIRs) of chikungunya (CHIKV), dengue (DENV), and Zika (ZIKV) viruses in Ae. albopictus and weather data (mean monthly relative humidity [RH] and temperature, and monthly total rainfall [RF]) in Cebu city during the dry (March -May 2021; March-May 2022) and wet (June-November 2021; February and June 2022) seasons *Statistically significant